% 周期问题的函数
function[c,f,s]=pdefun(x,t,u,ux)

c = [1;1];
% d1=0.21，d2=0.24, theta=1.7+0.12*abs(cos(pi * t)), b=0.4*abs(cos(pi*t)), c=0.1*abs(cos(pi*t)), m=0.884*abs(cos(pi*t)), r=2+0.2*abs(cos(pi * t)), f=0.14*abs(cos(pi*t))
%ρ(t)=1, 
%p = 1
%ρ(t)=1+0.4*abs(sin(pi*t)

p = 1+0.7*abs(sin(pi*t))

y(1) = (1.7+0.12*abs(cos(pi * t))) * u(1) - (0.4*abs(cos(pi*t))) *(u(1))^2 - (0.1*abs(cos(pi*t))) * u(1) * u(2) / (0.884*abs(cos(pi*t))+u(1));

y(2) = (2+0.2*abs(cos(pi * t))) * u(2) - (0.14*abs(cos(pi*t))) * (u(2))^2 / (0.884*abs(cos(pi*t)) + u(1));


F=[y(1);y(2)];

s=F.*[1;1];

f = [(0.21*ux(1))/(p^2);(0.24*ux(2))/(p^2)];
